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Abstract 

It is shown, concerning equivalent classes, that on a one-dimensional lat- 
tice with nearest neighbor interaction, there are only four independent 
models possessing double-shocks. Evolution of the width of the double- 
shocks in different models is investigated. Double-shocks may vanish, and 
the final state is a state with no shock. There is a model for which at large 
times the average width of double-shocks will become smaller. Although 
there may exist stationary single-shocks in nearest neighbor reaction dif- 
fusion models, it is seen that in none of these models, there exist any 
stationary double-shocks. Models admitting multi-shocks are classified, 
and the large time behavior of multi-shock solutions is also investigated. 
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1 Introduction 



Recently, shocks in one-dimensional reaction-diffusion models have absorbed 
much interest [1-18]. There are some exact results on shocks in one-dimensional 
reaction-diffusion models as well as simulations, numeric results [6] and also 
mean field results [2] . Formation of localized shocks in one-dimensional driven 
diffusive systems with spatially homogeneous creation and annihilation of parti- 
cles has been studied in [12]. Recently, in [4], the families of models with travel- 
ling wave solutions on a finite lattice have been presented. These models are the 
Asymmetric Simple Exclusion Process (ASEP), the Branching- Coalescing Ran- 
dom Walk (BCRW) and the Asymmetric Kawasaki-Glauber process (AKGP). 
In all of these cases the time evolution of the shock measure is equivalent to 
that of a random walker on a lattice with L sites with homogeneous hopping 
rates in the bulk and special reflection rates at the boundary [4]. Shocks have 
been studied at both the macroscopic and the microscopic levels and there are 
some efforts on addressing the question that how macroscopic shocks originate 
from the microscopic dynamics [7]. Hydrodynamic limits are also investigated. 

Among the important aspects of reaction-diffusion systems, is the phase 
structure of the system. The static phase structure concerns with the time- 
independent profiles of the system, while the dynamical phase structure concerns 
with the evolution of the system, specially its relaxation behavior. In [19- 
22], the phase structure of some classes of single- or multiple-species reaction- 
diffusion systems have been investigated. These investigations were based on the 
one-point functions of the systems. In a recent article both stationary and also 
dynamical single-shocks on a one-dimensional lattice have been investigated [18]. 
It was done for both an infinite lattice and a finite lattice with boundaries. Static 
and dynamical phase transitions of these models have been studied. It was seen 
that ASEP has no dynamical phase transition, but both BCRW and AKGP 
have three phases, and the system may show dynamical phase transitions [18]. 

The question addressed in this article is that, on a one-dimensional lattice 
with nearest neighbor interaction, which models possess double shock and also 
multi shock solutions. By double-shock it is meant an uncorrelated state where 
the occupation probability has two jumps. All the models have nearest neighbor 
interactions and are on a one-dimensional lattice. It is shown that, concerning 
equivalent classes, there are only four independent models possessing double- 
shocks. For two models, double shock disappears and the final state is a linear 
combination of Bernoulli measures. There is a model for which At large times 
the average width of double shock becomes small, stationary state is a state 
which does not evolve. It can be easily seen that there may exist stationary single 
shocks in nearest-neighbor reaction diffusion models (BCRW, and AKGP), in 
other words there are single shock states without any evolution. But in none of 
these models, there is no stationary double shock. Combining single shocks one 
may construct multi shocks. There are multi shocks of the type (0, p, 0, p, ■ ■ ■ ) 
and (0, 1, 0, 1, • • • ). At large times the final state is a linear combination of single 
shocks, or a state with no shock. 
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2 Notation 



Consider a one-dimensional lattice, each point of which either is empty or con- 
tains one particle. Let the lattice have L sites. An empty state is denoted by 
|0) and an occupied state is denoted by |1). 

10} := (T) , |1) := (f) • (1) 



If the probability that the site i is occupied is pi then the state of that is 
- / 
Pi 



represent by I 1 " <" I . The St a te of „,e sy ,e m . cH_ d by . ^ 



(2) 



where V is a 2-dimensional vector space. All the elements of the vector |P) are 
nonnegative, and 

(S|P> = 1. (3) 

Here (S| is the tensor-product of L covectors (s\, where (s\ is a covector the 
components of which (s Q 's) are all equal to one. The evolution of the state of 
the system is given by 

|P)=W|P), (4) 

where the Hamiltonian H is stochastic, by which it is meant that its nondiagonal 
elements are nonnegative and 

(S| H = 0. (5) 
The interaction is nearest-neighbor, if the Hamiltonian is of the form 

L-l 

H=J2 H i,i+i, (6) 

i=l 

where 

Hi i+ i := 1 <8> • • • <g> 1 ®H <g> 1 ® • • • <g) 1 . (7) 
v ' ^ ' 

j-l L-l-j 

Nondiagonal elements of shown by u>ij, are reaction rates, hence nonnegative, 
and its diagonal elements are nonpositive. u>ij is the rate for changes of the 
configuration of a pair of neighboring sites from the initial state j to the final 
state i. We take the state 1 00) as the state 1, |01) as 2, 1 10} as 3 and finally |11) 
as the fourth state. So, e.g. W23 is the rate for change of configuration 1 10) to 
1 01 ) , which is the hoping rate to the right. 

Any configuration of the system may be represented by the vector \E a ). So 
the system is spanned by 2 L vectors, \E a ) (a = 1,2, ■• - 2 L ), and any physical 
state is a linear combination of these vectors 

\F)=^V a \E a ), where = 1. (8) 

a— 1 a— 1 
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P a s are nonnegative real numbers. V a is the probability of finding the system 
in the configuration a. 

It is said that the state of the system is a single-shock at the site k if there 
is a jump in the density at the site k and the state of the system is represented 
by a tensor product of the states at each site as 



e k ) =u m ®v® {L ~ k \ (9) 



where 



It is seen that 



-i 1 ;") H 1 ;*)- (10) 

(S|e fc ) = 1. (11) 

\e k ) represents a state for which the occupation probability for the first k sites 
is pi, and the occupation probability for the next L — k sites is p?. The set 
\ek),k = 0, 1, • • • L is not a complete set, but linearly independent. 

There are three families of stochastic one-dimensional non-equilibrium lattice 
models, (ASEP,BCRW,AKGP), for which if the initial state of these models is 
a linear superposition of shock states, at the later times the state of the system 
|P) remains a linear combination of shock states. For these models 

H\e k ) = d\e k -i) + d'\e k +i) - (d + d')\e k ). (12) 

where d and d' are some parameters depending on the reaction rates in the bulk, 
and the densities p\ and p2- So the span of |efc)'s is an invariant subspace of W, 
the Hamiltonian of the above mentioned models. It should be noted that the 
number of \e k ) : s are L + 1, and an arbitrary physical state is not necessarily 
expressible in terms of |efc)'s. 

Let's assume that the initial state of the system is a linear combination of 
shock states 

L 

|P)(0) = $>fc(0)|e*>. (13) 

fe=0 

p k s, are not necessarily nonnegative, and so any of them may be greater than 
one. For such an initial state, the system remains in the sub-space spanned by 
shock measures. 

L 



Using (11), it is seen that 



|P>(t) = 5>*(*)|e*>. (14) 

k=0 



E»w = L ( 15 ) 



fe=0 

The three models are classified as following [4] 
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• 1. ASEP - The only non- vanishing rates in the bulk are the rates of diffu- 
sion to the right cj 23 and diffusion to the left w 32 . In this case the densities 
can take any value between and 1 (pi, p 2 ^ 0, 1). d, and d! are 

, pi(l-pi). . 

d = (w 23 - ^32) 

P2 - Pi 

a = (w 23 - w 32 ). (16) 

P2 - Pi 

It should be noted that the densities p\ , and p 2 are also related through 

P2(l- Pi) W 23 



Pl(l~P2) ^32' 

So 

Pi. . 7/ P2 . 



(17) 



d=^LJ 23 , d'^^32- (18) 

P2 Pi 

• 2. BCRW - The non- vanishing rates are coalescence (W34, and W24), Branch- 
ing (uj 42, and W43) and diffusion to the left and right (0)32, and ^23)- The 
density pi can take any value between and 1, but P2 should be zero. 
These parameters are related through 

^23 _ ^24 + ^34 _ 1 - Pi ^ 1Q ^ 
W43 ^42 + W43 p\ 

The parameters d, and d! are 

d = (1 - pi)w 32 + P1W34, 

<f = ^. (20) 
Pi 

if W32 = UJ34 = ^43 = ^23 = 0, and cj 24 /cj 42 = (1 — p)/p then = d' = 0, 
and the model admit stationary single shock. 

• 3. AKGP - The non-vanishing rates are Death (u>i 2 , and W13) and Branch- 
ing to the left and right (W42, and W43), and also diffusion to the left U32. 
pi should be equal to one, and p 2 should be zero. The hoping parameters 
are d = W13, d' = 0J43. 

3 Double shocks 

The state of a double shock may be defined through 

|e m>fc > = u® m ® v® k ® w® {L - k - m \ m + k<L, (21) 

where 

»-C;r) -C;r) -=-C;")- <»> 
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\&m,k) represents a state for which the occupation probability for the first m sites 
is pi, the occupation probability for the next k sites is p2, and the occupation 
probability for remaining sites is p 3 . We call such state a double shock, with 
the first shock at the site m, and the other one at the site m + k. k is the width 
of double-shock, and pi £ [0, 1]. To have a double shock pi should be different 
from p2, and P2 also should be different from p 3 . We search for Hamiltonians 
for which the span of |e m fc)'s is an invariant subspace of H, 

7~t\em,k) — rfi|e m -i,fe+i) + d[\e m+ i.k-i) + d 2 \e m .k-i) + d' 2 \e m ,k+i) 

-{d l + d' 1 +d 2 + d' 2 )\e m . k ), fc>2, (23) 

where rfj's (c^'s) are parameters depending on the reaction rates and may be 
considered as the rates of jump of the shock to the left (right). W|e m ,i) will be 
discussed later. 

As it is seen for single-shocks, one should study cases with different values of 
p separately. One may divide the region of values for p to p = 0, < p < 1, and 
p = 1. From now on the cases p = 0, and p = 1 will be explicitly stated, and 
whenever we write p, it is meant that p ^ 0, 1. To have a double-shock there may 
be different combinations of densities. There are different models, which may 
transform to each other through particle-hole, or right-left interchange. We call 
these models equivalent models. As an example the model admitting the double- 
shock (pi,P2,P3) = (0,p, 1) is related to the model admitting the double-shock 
(l,p, 0) through right-left interchange, and also related to the model admitting 
the double shock (1, 1 — p, 0) through particle-hole interchange. Let's consider 
the double shock (0, l,p). A necessary condition for the Hamiltonian for which 
the span of double-shock measures be an invariant subspace of 7i, is that the 
span of each of single-shock measures (0, 1) and (l,p) are separately invariant 
subspace of H. The single-shocks (0, 1) form an invariant subspace for the 
hamiltonian in AKGP. The only interactions which may have nonzero rates are 

%A -► (00, AA), A% -► (00, AA, %A). (24) 

As far as we consider the single shock (0, 1), there is no extra constraint on the 
nonzero reaction rates. The single-shocks (l,p) form an invariant subspace for 
the hamiltonian in BCRW, with the following interactions 

00 - (U, M), U - (00, A0), A0 -» (00, U), (25) 

whose reaction rates should satisfy 

V21 + ^31 _ ^23 _ P 
^12 + 0Ji 3 W13 1 - p 

The space of parameters of the model, for a double-shock (0, 1, p), is the overlap 
of the space of parameters of the AKGP and BCRW. Gathering all these together 
it is seen that all the reaction rates should be zero. So there is no reaction 
diffusion model with nearest neighbor interaction for which the double shocks 
(0, 1, p) form an invariant subspace. 

It can be easily shown that, concerning equivalent classes, there are only 
four independent cases. 
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• 1- (P1,P2,P3), 

Among the models possessing shock solution, (and for pi ^ 0,1), ASEP is 
the only model for which double shocks forms an invariant subspace. The 
only nonvanishing rates are U23 and W32, and they should satisfy 

^23 P2O--P1) p3(l~P2) 



^32 - Pi) /9 2 (l-/9 3 ) 

di's and rf-'s are 



(27) 



dl El= d2 El =(j23 ^ = ^ =W3a (28) 

Pi P2 P2 P3 

This model has been studied in [18]. UI23 and w 32 are positive nonzero 
rates. So dj's and d'^s are also nonzero. To have stationary double shock 
W|e TOi fe) = leading to dj = d\ = 0, which is unacceptable. So It is not 
possible to have stationary double shock in ASEP. 

2. (0,p,0), (or (p,0,p)) 

The necessary condition for a model possessing double-shocks (0, p, 0) (or 
(p, 0, p)) is that this model possesses both single shocks (0,p), and (p, 0). 
Nonvanishing rates for such a model are W23, ^24, ^32, W34, LO42 and W43. 
These rates should satisfy 

U>24 + W34 W32 W23 1 — P 



LO42 + W43 CJ42 ^43 



(29) 



cL 's and ci'-'s are 



dl = —, d[ = (1 - p)w 2 3 + P^24 
P 

d' 2 = —j d 2 = (1 - p)^32 + P^34- (30) 
P 

The Hamiltonian with the above mentioned reaction rates also possesses 
the double-shock (p, 0, p). The only difference is that the rate of jump to 
the left (and right) of the first double-shock is the rate of jump to the right 
(and left) for the second one. To have stationary double shocks d^s should 
be zero, which needs all the rates to be zero. So, there is no stationary 
double shock in this model. 

3. (0,1,0), 

Nonvanishing rates are W13, Wi 2 , U42, W43. This model is an asymmetric 
generalization of zero temperature Glauber model, di's and d'^s are 

di = LO42 d[ = U12 

d 2 = W13, d' 2 = W43. (31) 

To have stationary double shock di's and d-'s should be zero, which leads 
to vanishing all the reaction rates. So, this model does not have any 
stationary double shock either. 
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• 4. (0,p, 1). 

The only nonvanishing rate is u>23- di's and c^'s are 

di = 0, 4 = (1 - p)w 2 3 
d 2 -/9w 23 , 4-0. (32) 
This model does not have any stationary double shock either. 
If the initial state is a linear combination of double shocks, then 

oo oo 

l P )W= £ X>m*(*)|e ro *>, (33) 

m— — oo /c=l 

where p m ,k is the contribution of the double-shock, (mk) in the state of the 
system. Using (4), (33), and also the linear independency of |e mi fc)'s, one can 
obtain the evolution equation for p m k&- It is a difficult task to solve difference 
equations of this type. p m ^ has two indices, m representing the position of the 
first shock, and k the width of the double shock. We may forget about the 
position of the first shock and ask only about the width of double shock. Then 
one may encounter with a difference equation which can be solved more easier. 

Let's consider the general case where the span of |e m fe)'s is an invariant 
subspace of H 

m' ,k' 

If the Hamiltonian has the property that J2 m ' H™k i s independent of m, then 
one may define a new Hamiltonian TL through 

n k k--=J2 n ^ k ' (35) 

m' 

It can be easily shown that H. is stochastic, it is meant that 

H% > 0, for k' ^ k 

£ fe < Hi - 0. (36) 

Then one may forget about m, position of the first shock, and only ask about 
the contribution of double shocks with the width k. It is obvious that some 
part of information about the position of the first shock will be lost. Now one 
may define \fk) as the state of a double shock with the width k. Identifying all 
|e m ,fe) with the same m to each other in the state (33), one may define another 

state |P) where the information of the position of the first shock has being lost 

oo 

ip>(t) = X>(*)i/fc>- ( 3? ) 

k=l 
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Here qk(t) is defined through 

qk 



OO 

m=— oo 



(38) 



and it is the contribution of all double shocks with the width k. Then instead 
of (23), one may obtain 

H\f k ) = D\f k+1 )+D'\f k ^)-(D + D')\f k ), k>2. (39) 

where 

D:=di+d' 2) D':=d' 1 +d 2 . (40) 

3.1 Double-shocks (0, p, 0) and (0, 1,0) on a periodic lattice 

Let's consider a lattice with L sites and with periodic boundary conditions. 
Then, the only double shocks which could exist, are (0, 1, 0) and (0, p, 0). Let's 
sum up the contributions of all double shocks with the same width. The position 
of double shocks will be again washed out. Then one should work with \f k ), 
which stands for the state of a double-shock with the width k. /o), and 
are Bernoulli measures corresponding to an empty lattice and a full lattice, 
respectively. It can be easily shown that 



n\f ) = o, 

H\f k ) = D\f k+1 ) + D'l/fc-i) ~(D + D')\f k ) k 0, L, 

n\f L ) = o, (4i) 

Here D stands for the rate of increasing the width of double-shock, and D 1 
stands for the rate of decreasing its width. One can map this model to a model 
with one particle on lattice with boundaries at k = 0, and k — L. This particle 
hops to the right and left with the rates D and D', and there are traps at the 
boundaries. The system has only two stationary state, |/o), and //,), means 
that at large times there is no shock, and the final state is a linear combinations 
of the Bernoulli measures. 

W) = Qo\fo)+qL\fL). (42) 
If the initial state is a linear combination of |/fe)'s then 

F>(*) = X>(*)l/fc}- (43) 

Using (40), one arrives at 
q = D'qx, 

q\=D'q 2 -{D + D')q 1 , 
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q k = D'q k+1 + Dq k - X — (D + D')q k k ^ 0, 1, L - 1, L, 
Ql-i = DqL-2 — {D + D')qL-i, 

4l = DqL-i. (44) 
qk{t)'s in the bulk, (fc ^ 0, L), can be obtained. They are 

»(*) - |(§) fc/2 e - (i?+D ' )t i:E^( o )(#) m/2 

s=l m— 1 

sin (^) s i n (!^) e ^co S (Wi). (45) 
Li ij 

One may integrate qi(t), and qL-i(t) to obtain g'o(t), and <7i(t), which are the 
only terms surviving at large times. There is also another way to obtain the qo, 
and <jx at infinitely large times. In fact, there are two constants of motion X\ 
and 1% ■ T\ is related to the conservation of probability 

L 

(S|P) = 1 => Ji:=2« fc (t) = l, (46) 

fe=0 

and 

fc=0 

It should be noted that the system has two stationary states, so there are 
two right eigenvectors corresponding to zero eigenvalue for the Hamiltonian Tt. 
Therefore there are also two left eigenvectors corresponding to zero eigenvalue 
for H. These are 



(S| = (l 1 1 ••• 1). (48) 
<S'|=(l ^ (^) 2 (^) 3 ••• (49) 

The second constant of motion can be obtained using (§'|P)(£). As long as 
D ^ D', the constants of motion X\ and 2" 2 are two independent quantities. 
For D = D', X\ and 1 2 arc the same. But as the stationary state has two- 
fold degeneracy, there should exist another constant of motion. The second 
independent constant of motion is I' 2 ■= J2k=o ^ Qk(t) = (k)- So, for D = D', 
the average width of the shock, (k), is a constant of motion. One should expect 
this, because D and D' are the rates for increasing the width of the double 
shock and decreasing it respectively. 

For the double-shock (0, p, 0), D' / D = 1 — p < 1. So, the constant of motions 
are X\ and Z2 ■ X\ is the summation of probabilities for finding a double shock 
with any width, so it should be equal to one. The second constant of motion 
also has a physical interpretation. The rate for changing any configuration 
of a pair of neighboring sites to the state 1 00) is zero. So the probability for 
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finding a completely empty lattice does not change with time. I2 is exactly the 
probability of finding an empty lattice in the initial state 



(50) 



k=0 



Using constant of motions, for D 7^ D' , at infinitely large times, we have 

qo + qL = l, 



?o + (§)V = X>(0)(§) fe - 



Solving these equations one obtains 

L 

qo( 



(51) 



)(oo)= [£<fe(0)(^) fc -(^ 



fe=0 



1 - ( — ) L 



gi (oo)= [l-£ft(0)(^) 



fc=0 



(52) 



As it is seen the contribution of |/o) and in the final state depends on both 
reaction rates and initial conditions. 

The Hamiltonian for the model possessing the double-shock (010), with D = 
D' is the Hamiltonian for zero temperature Glauber model. This model have 
been studied in [19,23,24]. The average density at each site (iii)(t) at the time 
t and also all the correlation functions at large times for an infinite lattice have 
been calculated in [23] . Static and dynamical phase transitions of this model 
have been also studied in [19]. Here, D is not necessarily equal to D' . For 
D' > D, and large L, one arrives at 



D' 



Qo 



( oo) = l-E%(0)(^) 



k-L 



k=0 



<7l(oo) 



9o(oo), 



(53) 



If initially only double shocks with finite widths have contributions, in the ther- 
modynamic limit (L — ► 00) the system will finally fall in the state /o, but if 
D' < D it can be seen that both stationary states have contributions in the 
final state. 

For the case D = D' one obtains 



1 L 1 
3o(oo) = l--Y J kq k {Q) = l--{k), 

k=0 

1 L 1 

9l(oo) = yE fc 9fe (°) = 7W' 



(54) 



k=0 
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which means that at large times the system is fully occupied or empty. The 
probability of finding a fully occupied lattice at large times is equal to the 
initial average width of the double-shock divided by the size of the lattice. 



3.2 Double-shock (0, p, 1) 

Let's consider the double-shock (0,p, 1) on an infinite lattice. The only non- 
vanishing rate is u>23, which can be set equal to 1, by a suitable redefinition of 
time. Direct calculation gives 



H\e m ,i) = 0, 

H\e m ,k) = (1 - p)|e TO +i,fc-i) + P \e m ,k-i) - \e m ,k), k ^ 1. (55) 

It is seen that there is no probability for width increase. If there is initially a 
shock |e m fe), at later times its width becomes smaller, and at large times there 
are only double shocks with the width 1. Starting with a linear combination of 
the shocks, the evolution equation for p m fc's, can be obtained to be 

Pm,l = (1 - p)p m -l,2 + PVm.2, 

Pm,k = (1 - P)Pm-l,k+l + PPm,k+\ ~ Pm,k fc ^ 1. (56) 

Defining q k , through (38), one arrives at 

<7i = Q2, 

Qk = Qk+i - qk k^l. (57) 

If initially the state of the system is a double shock, e.g. \cmk)- Then, it 
is obvious that at later times there are only double shocks with the position 
of the first shock in the range M < m < M + K — 1, and with the width 
1 < k < M + K — m. Let's assume the initial state is 

L 

|P>=I>(0)|/k>, (58) 

fc=0 

where \fk) is the state of double shocks with the width k. Then 

<?o = 0, 
qi = <?2, 

4k = qk+i - qk, 2 < k < L - 1, 

qh = ~qL, 

4=0, L + l<k. (59) 

The above equations show that at large times there are only contributions of 
the double shocks with the width 1. These set of equation can be solved leading 
to 

L — k n 

q k (t) = J2 1k+n(0)— e-* 2 < k < L. (60) 

n=0 



11 



This together with q\ = q 2 can be used to obtain qi(t). 

ft ,/ n 



Ql(t) = «l(0) + V q n +2(0) / -r e"*' 

L L— 2 n m 

^^"EE^O^e-'. (61) 

n— 1 n— m=0 

As it is expected at large times all the double shocks changes to the double-shock 
with the width 1. 

Let's study the distribution of these double shocks at large times. Using 
(55), and defining A m . k := cxp(tH)\e mk ) , it is seen that 

^4^— + A m<k =(1 - p)A m+ i <k _ 1 + pA m ^ k _i k^l, 

A mA =\e kl ). (62) 
At large times this equation recasts to 

A TOife (oo) = (1 -p)A m+ i )fc _i(oo) +pA TOife _i(oo) fc^l, (63) 
whose solution is obtained to be 

A m . fe M = lim (e tH \e mk )) = ]T ( fc 7 ^ (1 - p)'^ -1- '' |e ro+j ,i). (64) 

j=0 \ J / 

So, at large times the state of the system is a linear combination of double 
shocks with the width 1. The distribution of the position of these double shocks 
is a binomial distribution. Let's consider the initial state to be a double shock 
with the width k, |eo,k), then the average position of the first shock at large 
times is 

(j) = (fc-l)(l-p), (65) 



and the width of the binomial distribution is \J p(l — p)(k — 1). 

4 Multi shocks 

Combining single shocks one may construct multi shocks. The only models with 
multi shocks are as following 

• L (pi,P2,P3,---) 

The span of multi shocks (pi, p2, P3, • • • ) is an invariant subspacc of Hamil- 



tonian of ASEP provided the densities satisfy 

P»+i(l - Pi) _ ^23 
Pi(l-p i+1 ) LU 32 ' 



(66) 
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It can be easily seen that the rate of hoping of the zth shock to the left, 
di, and also the rate of hoping of the ith shock to the right, d-, is given by 



di =U>23 



Pi 




Pi+l 
Pi+l 

Pi 



(67) 



This result are obtained in [3]. 

• 2.(0,p,0,p, •••) and (0, 1, 0, 1, • • • ) 

The model admitting multi shock of the type (0, p, 0, p, ■ ■ ■ ) are that ad- 
mitting double shocks (0,p, 0), or (p,0,p). The model Possessing multi 
shocks of the type (0, 1,0, 1, • • •) is an asymmetric generalization of the 
zero temperature Glauber model. In both of these multi shocks there are 
edges at shock points. The edges are destroyed two by two. Let's consider 
a multi shock of order N with the first shock at the site m. It is seen that 
the action of Hamiltonian on such state is 

H\em,ki,— ,fciv-i} = rfi|e m -i,fci+i,--- : fe JV _ 1 ) + dilem+i^-i,... .fcjv-i) 



If any of the fc^s in the left hand side is equal to one, at the right hand 
side there will be a multi shock of order N — 2. So there is a finite 
probability that the system falls in a state with lower shocks, and there is 
no probability for increasing the number of shocks. In fact if there exists 
a state for which any state can transform directly or even indirectly to 
it, and that state has no evolution, then that state is the final stationary 
state. Let's consider a periodic lattice. The number of shocks, N, should 
be even. So at large times the state of system is a state with no shock. 
For the models on an infinite lattice number of shocks, N, may be even 
or odd. Then For odd N the final state is a linear combination of single 
shocks. 

5 Summary 

There are three types of models with travelling wave solutions on a onc-dimcnsional 
lattice. These are classified in [4]. It is seen that there are four type of models 
admitting double shocks. Double shocks and the models admitting these double 
shocks are as following 

• {pi,P2,Pz)- Nonvanishing rates are W23, and W32. 

• (0, p, 0), (and also (p, 0, p)). Nonvanishing rates are W23, ^24, ^32,^34, ^42 



+d2|e m ,fe 1 -i,... ,fe 



: JV-1 



) +4l e m,/ Sl +l,-,fe N _ 1 ) H 



+d 2 \ > + 41 ). (68) 



and W43. 
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• (0, 1, 0). Nonvanishing rates are ^13,^12,^42,1^43. 

• (0, p, 1). The only nonvanishing rate is is cj 2 3- 

There are three type of models admitting multi shocks. The multi shocks 
are of the type 

• (pi, Pi, Pi, ■ ■ ■ )• Nonvanishing rates are o>23, and a> 32 . 

• (0, p, 0, p, ■ ■ ■ ). Nonvanishing rates are W23, ^24, ^32, ^34, ^42 and W43. 

• (0, 1,0, 1, • • ■). Nonvanishing rates are W13, W12, w 42 , W43. 
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